Testing intermediate-age stellar evolution models with VLT 
photometry of LMC clusters. II. Analysis with the Yale models 
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ABSTRACT 

We present an analysis, using Yale stellar evolution models, of the color-magnitude diagrams 
(CMDs) of three intermediate-age LMC clusters, namely NGC 2173, SL 556 and NGC 2155, 
obtained with the VLT. The main goal of our project is to investigate the amount of convec- 
tive core overshoot necessary to reproduce the CMDs of relatively metal-poor, intermediate age 
stellar populations, to check whether the extrapolation that is usually made from solar metallic- 
ity is valid. In the process, we obtained values for the binary fraction of each cluster, together 
with refined age estimates. Our method involved the comparison of the observed CMDs with 
synthetic CMDs computed using various values of the overshoot parameter and binary fraction. 
We conclude that a moderate amount of overshoot and some fraction of binary stars are essen- 
tial for reproducing the observed shapes around the turnoff in the CMD's of all three clusters: 
unresolved binary stars fill in the expected core contraction gap, and make a unique sequence 
near the gap, which cannot be reproduced by single stars alone, even with a larger amount of 
overshoot. We utilize ratios of the number of stars in different areas around the core contraction 
gap to constrain the binary fraction, which is around 10-20% (for primary-to-secondary mass 
ratio > 0.7) in all three clusters. Even if binary stars contaminate the core contraction gap, 
it is shown that the overshoot parameter can be inferred from the color dispersion of the stars 
around the contraction gap, regardless of the assumed binary fraction. From our overall analysis 
such as, shape of isochrones, star counts, color distribution, and synthetic CMD comparisons, 
we conclude that overshoot ~ 20% of the local pressure scale height best reproduces the CMD 
properties of all three clusters. The best age estimates are 1.5, 2.1 and 2.9 Gyr for NGC 2173, 
SL 556 and NGC 2155, respectively. 



Subject headings: 
evolution 



color-magnitude diagrams — galaxies: star clusters — Magellanic Clouds — stars: 



^Department of Astronomy, Yale University, P.O. Box 
208101, New Haven, CT 06520-8101; jhwoo@astro.yale.eu, 
demarque@astro.yale.edu 

^ Andes Prize Fellow, Universidad de Chile and Yale Uni- 
versity 

^Currently: Ramon y Cajal Fellow. Instituto de As- 
trofi'sica de Canarias, 38200 La Laguna, Tenerife, Canary 
Islands, Spain; carme@iac.es 

^University of Oxford, Astrophysics, Keble Road, Ox- 
ford OXl 3RH, UK; yi@astro.ox.ac.uk 

® European Southern Observatory, Karl-Schwarzschild- 
Strasse 2, D-85748 Garching bei Miinchen, Germany; 
mzoccali@eso.org 

Based on observations collected at the European 
Southern Observatory, Chile (ESO N° 64.L-0385). 



1. Introduction 

Stellar evolution theory is routinely tested us- 
ing Galactic clusters, either globular (in general 
old, and both metal-poor and metal-rich) or open 
(intermediate-age or young, and generally more 
metal-rich than globular clusters). These clusters 
reflect the particular star formation and chemi- 
cal enrichment history of the Milky Way, and cer- 
tainly do not represent all possible stellar popula- 
tions. For example, intermediate-age, metal-poor 
populations, which are not well represented in the 
Galaxy, are conspicuous in the Magellanic Clouds. 
This is an important part of the parameter space, 
which is relevant, for example, to the study of 
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stellar populations of dwarf galaxies in the Local 
Group, or to that of galaxies at high redshifts, and 
therefore, in early evolutionary stages. 

There are uncertainties that are critical to 
modeling intermediate-mass stellar models. Most 
notable of them are the treatments of convec- 
tion (e.g., convective core overshoot) and mass 
loss. In this paper, we will address them, us- 
ing intermediate-age Magellanic Cloud clusters as 
comparison templates. 

Convective core overshoot (OS) is one of the 
physical parameters which must be taken into ac- 
count during the core hydrogen burning phase of 
the evolution of intermediate-mass stars. Stellar 
evolution theory predicts that the CMD of young 
and intermediate-age star clusters is greatly af- 
fected by the amount of OS. In the early days 
of stellar structure theory, the convective core 
size was determined by the classic Schwarzschild 
(1906) criterion, ignoring the so-called OS, which 
is due to the inertial motion of materials beyond 
the formal core edge. Initially, the extent of 
core OS was thought to be negligible (Saslaw & 
Schwarzschild 1965). Subsequent theoretical stud- 
ies, however, have emphasized the importance of 
convective core OS in stellar evolutionary models 
(Shaviv & Salpeter 1973). 

Since then, many observational and theoretical 
studies have been devoted to investigating the ef- 
fects of core OS (Prather & Demarque 1974; Bres- 
san, Chiosi & Bertelli 1981; Bertelli, Bressan & 
Chiosi 1985; Stothcrs 1991; Mcynet, Mcrmilliod 
& Maeder 1993). It is now generally believed that 
a moderate amount of core OS is necessary in or- 
der to reproduce the shape of the core contrac- 
tion gap observed in CMDs of galactic open clus- 
ters (Maeder & Mermilliod 1981; Stothers & Chin 
1991; Carraro ct al. 1993; Daniel et al. 1994; 
Demarque, Sarajedini, & Guo 1994; Kozhurina- 
Platais et al. 1997; Rosvick & VandenBerg 1998), 
but the dependence of core OS on stellar mass 
is still poorly understood and the subject of on- 
going investigations (e.g. Ribas, Jordi, & Gimenez 
2000). 

The extent of OS is not currently derivable 
through direct theoretical calculations but can 
be determined empirically by comparing synthetic 
CMDs to observed ones, with OS in the models of- 
ten being parameterized as a fraction of the pres- 
sure scale height at the formal edge of the convec- 



tive core. In the CMDs of intermediate-age star 
clusters (2 < Age < 5Gyr), which have a dis- 
tinctive gap near the turnoff (TO) caused by the 
rapid contraction of a hydrogen-exhausted convec- 
tive core, the TO topology can be used to evalu- 
ate the extent of core OS (Demarque et al. 1994; 
Kozhurina-Platais et al. 1997). In the case of 
younger clusters, the importance of core OS can 
be estimated from the observed luminosity func- 
tion of main sequence (MS) stars (Chiosi et al. 
1989; Barmina, Girardi & Chiosi 2002) and from 
the ratio of main sequence to post-main sequence 
stars. Studies of detached eclipsing binaries can 
also provide empirical estimates of OS by direct 
comparison with stellar tracks in the log Tes - log 
g plane (Schroder, Pols, & Eggleton 1997; Guinan 
et al. 2000; Ribas et al. 2000). However, the num- 
ber of binary systems available for OS studies is 
very limited. 

Many galactic open cluster studies have sug- 
gested a moderate amount of OS, i.e. ~ 20% of 
the local pressure scale height at the edge of a con- 
vective core. However, this has been tested only 
for near solar metallicity clusters in the Galaxy. 
So far, the presence of OS in metal-poor clusters 
in the Magellanic Clouds has only been tested in 
a limited sample of young (age < 500 Myr) clus- 
ters (Lattanzio et al. 1991; Vallenari et al. 1991, 
1994; Stothers k Chin 1992; Brocato, Castellani 
& Piersimoni 1994; Chiosi et al. 1995; Testa et al. 
1999; Keller, Da Costa & Besseh 2001; Barmina et 
al. 2002). For Magellanic Clouds intermediate-age 
clusters, the necessary observations are more chal- 
lenging due to the fainter magnitude of the TO, 
and HST or 8-m class telescopes under excellent 
seeing conditions are required. 

The characteristics of mass loss during the RGB 
phase are still unclear in many aspects. Among 
the early papers on the subject, that of Reimers 
(1975) was most influential and his mass loss for- 
mula has been widely used to date: 

M = -4 X 10-"r/— (1) 
gR 

where L, g, and R are luminosity, gravity and 
radius, respectively. The mass loss efficiency pa- 
rameter T], which was later introduced, has been 
reported to vary from 0.25 to 2 - 3 for metal-rich 
stars (Dupree 1986; Kudritzki & Reimers 1978; 
Renzini 1981). In contrast, the estimated range 
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of 77 for nictal-poor stars appeared to be nar- 
row, as most studies on metal-poor stars sug- 
gest 77 = 0.3 - 0.7 (Aaronson k Mould 1982; 
Mould & Aaronson 1982; Renzini 1981; Lec ct al. 
and these values successfully reproduce the hor- 
izontal branch morphology of the old Galactic 
globular clusters (Yi ct al. 1999). While the ap- 
plication of Reimers' formula for low mass stars 
predicts that the total amount of mass loss gets 
smaller as mass increases, it is not clearly demon- 
strated how the amount of mass loss depends on 
the mass of the stars in intermediate-age clusters 
(Tripicco, Dorman & Bell 1993; Liebert, Saffer & 
Green 1994; Carraro et al. 1996). 

In this context, we started a project to in- 
vestigate several aspects in the evohition of 
intermediate-age, metal-poor stellar populations, 
in particular, the amount of OS and the mass 
loss characteristics. To reach this goal, we ob- 
served three intermediate-age clusters using the 
Very Large Telescope (VLT) under excellent see- 
ing conditions. The project and the observations 
are presented in Gallart et al. (2002; Paper I). In 
the current paper, we will compare the observa- 
tions with the recent set of Yonsei-Yale (F^) stellar 
evolutionary models, and with models specifically 
calculated for this project using the same main 
input physics (but differing in the OS parameter), 
to constrain the amount of OS and the mass loss. 

The three clusters under the study seem to 
present the signature of a core contraction gap in 
their TO region, which is predicted by the mod- 
els as well for this age range. A possible obstacle 
in studies of OS using the shape of the core con- 
traction gap is binary contamination. Castellani, 
DegTInnocenti & Marconi (1999) pointed out that 
the contraction gap enhanced by core OS can be 
filled in by unresolved binary stars. In fact, some 
intermediate-age clusters, which are expected to 
have a core contraction gap, show a continuous 
sequence rather than a distinctive gap (Carraro et 
al. 1994; Paper I). Although this could be consid- 
ered as a limiting factor for the topology method, 
we will show that it is possible, using synthetic 
CMDs, to put some constraints on the binary star 
population. It may be also expected that an in- 
crease of binary fraction in synthetic CMD simula- 
tions tends to decrease the value of OS required for 
a good fit. Testa et al. (1999) have claimed that 
a no-OS model with ~ 30% of unresolved bina- 



ries allows a fair fit to the observed MS luminosity 
function of NGC 1866 in the LMC, while Barmina 
et al. (2002), using a revised analysis of the same 
, data, have concluded that the inclusion of unre- 
solved binaries would not alter the need for OS 
(see also for other clusters, Vallenari et al. 1992; 
Carraro et al. 1994) 

In this paper, we present a short description of 
the observational data (§ 2). Stellar models in- 
cluding OS treatment and synthetic CMD simula- 
tions are described in section 3. Detailed analysis 
on the steps of deriving ages, binary fraction, the 
amount of OS and mass loss is described in sec- 
tion 4, followed by a short discussion (§ 5) and 
conclusions (§6). 

2. THE DATA 

The three LMC clusters NGC 2173, SL 556 and 
NGC 2155 were observed with FORSl at VLT- 
UTl, in service mode, through the filters V and 
R, under excellent seeing conditions. PSF fit- 
ting photometry was obtained with the suite of 
programs DAOPHOTII/ALLSTAR/ALLFRAME 
(Stetson 1987, 1994), and calibrated using data 
of the same clusters obtained under photomet- 
ric conditions using the MOSAIC camera at the 
CTI0-4m telescope. We estimate an accuracy of 
the zero points of the photometric transformation 
of ± 0.005. An assessment of the completeness 
and error characteristics of the data was performed 
through artificial star tests. Finally, we carefully 
subtracted the LMC field stars in the cluster field 
in a statistical fashion. All these aspects of the 
data treatment are discussed in detail in Paper I. 
In the following sections we will just use the photo- 
metrically calibrated, field-subtracted CMDs, and 
the results of the artificial star tests for complete- 
ness and errors. 

3. MODELS 

3.1. Stellar tracks & isochrones 

We have constructed stellar tracks of metallicity 
Z = 0.004 and 0.007 with the Yale stellar evolu- 
tion code YREC (Gucnther et al. 1992), ignoring 
the effects of rotation. All input physics and model 
construction assumptions are the same as in the 
y2 Isochrones (Yi et al. 2001) except for the OS 
treatment. Three sets of tracks have been gener- 
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atcd depending on the amount of core OS. The OS 
parameter is defined as a fraction of the pressure 
scale height at the edge of a convective core and 
prescribes the radial size of the OS region beyond 
the classical core convection boundary. We used 
10, 20 and 30% of the local pressure scale height 
as the amount of OS, corresponding to "OS pa- 
rameter" values of 0.1, 0.2 and 0.3, respectively. 

Prom these tracks, we have constructed isochrones 
between 1 and 3Gyr, with an interval in age of 
0.1 Gyr. Color transformations are based on Leje- 
une, Cuisinier & Buser (1998), following the 
scheme. It is generally found that a convective 
core develops when the stellar mass is larger than 
~ 1.1 Mq, a result known to depend on chemical 
composition. For Y"^ isochrones, Yi et al. (2001) 
have adopted OS = 0.2 for young isochrones (< 
3 Gyr) and OS = 0.0 for old ones (> 3 Gyr) based 
on the observational studies of approximately so- 
lar abundance populations. Their isochrones have 
a 1 Gyr age grid for ages greater than 1 Gyr and 
do not provide the isochrones of ages greater than 
2 Gyr but smaller than 3 Gyr. This is inconve- 
nient to our study because our sample clusters 
appear to be in this age range. Besides, the 
Isochrones group does not recommend any sim- 
plistic isochrone interpolation for the age range of 
2 to 3 Gyr because it is in this age range that the 
transition between convective core and radiative 
core in the MS TO stars occurs (see also Woo & 
Demarque 2001). For this reason, we have com- 
puted a new Y'^ Isochrones in which the OS pa- 
rameter has been kept constant for ages between 
1 and 3 Gyr. Figure 1 illustrates the isochrones of 
two different OS parameters. Models with larger 
OS show brighter TOs and subgiant branches for 
each given age, implying that an age estimate for 
an observed CMD would be larger when using 
isochrones with larger OS. 

3.2. Synthetic CMD code 

Stellar mass is the fundamental parameter that 
determines the evolutionary path of each star in 
the CMD. By randomly distributing stellar mass 
with an initial mass function (IMF), and tracing 
each star's evolutionary stage at a given time, one 
can construct synthetic CMDs for stellar systems 
with single age and metallicity. These simple stel- 
lar populations are relatively easy to model due to 
the simplifying assumption of instantaneous star 



formation rate and homogeneous chemical com- 
position. 

Several kinds of uncertainties still remain in 
computing synthetic CMDs for simple stellar pop- 
ulations. The first kind, which affects each star's 
location on the CMD, comes from uncertainties in 
stellar evolutionary theory and stellar atmosphere 
models. Secondly, uncertainties in the number of 
stars in each location on the CMD come from the 
combinations of uncertainties in the evolutionary 
time scale of each stage, the slope of the IMF, and 
stochastic effects caused by small number statis- 
tics. Many CMD studies have aimed at testing 
physical parameters in stellar evolutionary theory 
by comparing theoretical shape and star counts 
with the observations while minimizing uncertain- 
ties from other sources. Depending on the nature 
of the targets and on the methods used in the anal- 
ysis, these uncertainties can be minimized with 
rigorous effort. For example, plausible stochastic 
effects in number ratio analysis can be avoided by 
repeating the simulations and using large enough 
areas in the CMD to count stars in observed and 
synthetic CMDs. 

A new CMD synthesis code has been created, 
including binary simulation, based on previous 
works (Woo 1998; Yi et al. 1999). For MS to 
RGB stars, Salpeter's IMF with x = 1.35 has been 
used. For synthetic horizontal-branch construc- 
tion, we have used the core helium-burning stellar 
models of Yi, Demarque & Kim (1997) with Lee, 
Demarque & Zinn (1990)'s prescription for mass 
distribution. Mass loss along the RGB is formally 
calculated with Reimers' (1975) formula, and sub- 
tracted at the RGB tip to derive the mass of core 
helium burning stars. However, for this particular 
study, the amount of mass loss is empirically de- 
termined to reproduce the observed luminosity of 
the red clump. 

Unresolved binary stars make a distinct se- 
quence in the observed CMDs. It is critical to in- 
clude binary stars for a detailed study of the CMD. 
We define the binary fraction (BF) as the number 
ratio of unresolved binary stars to the total num- 
ber of stars in the CMD. The binary stars in the 
CMD are simulated in the following way. Once 
stellar mass is distributed according to the stel- 
lar IMF, some stars are randomly selected to be 
given a secondary star. The mass of the secondary 
star Mse is assigned using the critical mass ratio. 
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Qc that sets the minimum mass ratio between sec- 
ondary and primary: 

M^e = {RAN * (1 - Qc) + qc) * Mpr (2) 

where Mpr is the primary star mass, and RAN is 
a random number with a flat distribution between 
and 1. Such a flat distribution for the mass 
function of the secondary stars and a typical value 
of ^0=0. 7 is adopted following Elson et al. (1998; 
see also Hurley & Tout 1999). Secondary stars 
with masses Mg^ < 0.7Mp,. would form a sequence 
indistinguishable from that of single stars. The BF 
that we will quote throughout the paper will be, 
therefore, the fraction of binary stars with Mge > 
0.7Mpr, and a larger total BF should be expected. 

Once the mass of the secondary star is deter- 
mined, then the magnitudes My(pr) and Mv'(se) 
of the primary and secondary, respectively, are 
combined to calculate the magnitude My(BI) for 
the unresolved binary system: 

Mv{BI) = -2.5*log [lO-0-i*Mv(pr)^-^^Q-OA*Myise)^ 

(3) 

3.3. Error &c completeness simulation 

The simulation of the observational errors (pho- 
tometric shifts and completeness) in the synthetic 
CMDs has been performed on a star-by-star basis, 
using an empirical approach without any assump- 
tion on the nature of the errors or their propaga- 
tion. The whole process has been performed in a 
way similar to that described in Aparicio & Gal- 
lart (1995) and Gallart et al. (1996a,b), and it 
is based on the results of the artiflcial star tests 
discussed in Paper 1. 

We have derived two tables from the artificial 
star tests. In the first one, the completeness in- 
formation for each V, R magnitude interval is 
recorded. In the second table, for each artificial 
star that was recovered under the same quality 
criteria as the cluster data, the input and out- 
put magnitudes are listed. From this informa- 
tion, the simulation process for photometric errors 
and completeness in the synthetic CMD is the fol- 
lowing: for each star in the synthetic CMD, with 
magnitudes Vs, Rg, we first use the completeness 
information in the corresponding magnitude inter- 



val to determine the probability of the star hav- 
ing been lost in the photometry, and keep it or 
remove it from the list of synthetic stars accord- 
ing to that probability. In a second step, for each 
synthetic star that has been kept, we select all ar- 
tificial stars within a given distance in color and 
magnitude from the synthetic star. One of these 
stars is then picked up at random, and its dV and 
5R, calculated as the differences between the re- 
covered and injected magnitudes for that star, are 
added to the magnitudes Vg, Rg in order to simu- 
late the observational errors. 

In Figure 2, synthetic CMDs before and after 
error simulations are compared as an example. In 
the left panel, single stars are placed on top of a 
given isochrone while the binary sequence has a 
dispersion according to the range of the mass ra- 
tio q ((/c < 9 < !)• The error simulated CMD 
in the right panel shows many scattered stars in 
the CMD due to the photometric errors especially 
at faint magnitude levels. The result of the com- 
pleteness simulation is that, globally, 10 to 15% of 
the total input stars are lost in the error simulated 
CMD. 

4. Analysis: overshoot and mass loss 

In this section, we compare the CMDs of the 
three LMC clusters, NGC 2173, SL 556 and 
NGC 2155, with isochrones and synthetic CMDs, 
in order to test the predictions of stellar evolution 
theory under different sets of assumed parameters 
regarding OS and mass loss. This will allow us 
not only to derive the cluster parameters (ages, 
metallicity, distance, and reddening) but also to 
constrain the input parameters in stellar model 
calculations (e.g., OS and mass loss). 

4.1. The amount of overshoot in low- 
metallicity, intermediate-age popula- 
tions 

Cluster CMDs provide a powerful tool for de- 
riving the OS parameter empirically. We have ap- 
plied the following strategy on our sample clus- 
ters. First, we used the color difference A(V-R) 
between the base of the RGB and the TO to con- 
strain the age range for each cluster under each OS 
assumption (Sarajedini & Demarque 1990). This 
has been complemented with visual isochrone fit- 
ting within that age range, to obtain a single best 
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fitting age for cacli OS parameter. Second, since 
the presence of unresolved binary stars substan- 
tially affects the morphology of the TO and core 
contracition gap in the CMD, we attempt to find 
the BF using the ratio of the number of stars in 
two boxes around and below the core contraction 
gap. The last step takes advantage of the differ- 
ent spread in color near the core contraction gap 
depending on the assumed OS, to single out the 
best OS, and consequently, the best age for the 
cluster. Binary stars have substantial effects on 
the morphology of the TO and the core contrac- 
tion gap, however, we will show that they do not 
greatly affect the color distribution of the stars, 
and therefore, the final conclusion on the best fit- 
ting OS parameter is almost independent of the 
assumed BF. These steps are explained in detail 
in the following three sections. 

4-1. 1. Isochrone solutions 

The major effect of core OS to stellar evolution 
is to supply hydrogen fuel from the OS-induced 
mixing region to the core, which lengthens the 
hydrogen core burning phase. This results in a 
brighter TO and subgiant branch since for a given 
age, more massive stars are going through the core 
contraction phase. Consequently, if the amount of 
OS assumed in the models is larger, an older age 
will be derived for a given observed CMD. There- 
fore, the OS parameter cannot be uniquely deter- 
mined without an accurate age constraint and vice 
versa. Due to this degeneracy, we determined the 
best age for each OS parameter in the following 
steps. 

As the TO gets redder with increasing age, the 
color difference, A(V-R) between the base of the 
RGB and the TO decreases. For a given age, 
a larger OS implies an increase of A(V-R) since 
more massive stars are located near the TO (Fig- 
ure 3). Thus, for each OS parameter, an age can 
be estimated from the observed A(V-R) value. 
However, due to the uncertainties in determining 
the locations of the base of RGB and the TO in 
the observed CMD, this step only gives a plausible 
range of age for each OS. A typical error, 0.02 in 
A(V-R) gives age uncertainty ~ 0.3 Gyr. 

Within the age range, we have performed 
isochrone fitting to find the best fitting age. In 
this step, two further parameters are required: 
distance modulus and interstellar reddening. The 



distance modulus of the LMC is usually taken to 
be 18.5 ± 0.1, which is consistent with the recent 
compilation of Benedict et al. (2002). However, 
there is still some disagreement on the LMC dis- 
tance modulus among experts in the field, and 
values ranging from 18.3 to 18.7 have been ad- 
vanced. In addition, the LMC has a large angular 
distribution of clusters in the sky, among which 
the distance modulus may differ. For example, 
if we consider that the line of sight depth of the 
LMC is of the order of the LMC bar size, 3 degree 
(~ 2.6 kpc), the distance modulus can easily differ 
by up to 0.1 magnitude. Currently, the distance 
to individual clusters is unknown. 

A mean reddening of E(B-V)~0.1 [E(V-R)= 
0.08] or larger is tipically found for the LMC (Har- 
ris, Zaritsky, & Thompson 1997, and references 
therein). Most of these studies, however, focus on 
the central part of the LMC, or on regions contain- 
ing young stars. Therefore, presumably, a larger 
amount of dust is present in these regions than 
in more external LMC regions, as is the case of 
the clusters in our study. Schlegel, Finkbeiner, 
& Davis (1998) calculate typical reddenings to- 
ward the LMC from the median dust emission in 
surrounding anuh, and find E(B-V)=0.075 (E(V- 
R)=0.06)), which should be close to the reddening 
found in the outer part of the LMC where internal 
extiction may be very low. No estimates of indi- 
vidual reddenings toward the clusters in our study 
are available. 

Given these uncertainties, we performed isochrone 
fitting with limited freedom of distance modulus 
and reddening. We did not use a fiducial line ob- 
tained from the CMD in the fitting process since 
a fiducial line is difficult to determine around the 
convective core contraction gap, and binary stars 
can influence the fiducial line position. Since field 
star subtraction is not perfect, we rather put more 
weight on the isochrone shape near the TO and 
upper MS, which are relatively well populated. 

With isochrone fitting, we single out the clus- 
ter age for each OS parameter, which is also com- 
plemented by comparing the observed CMD with 
the synthetic CMD simulated with each set of age 
and OS (Table 1). The true metallicities of our 
sample clusters are still poorly known but often 
suggested to be in the range Z=0.004 - 0.007 (Pa- 
per 1). We prefer to use a single value of metal- 
licity for all three clusters in order to investigate 
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systematic effects. Metallicity Z=0.007 seems to 
be inadequate as a representative value because 
the Z=0.007 isochrones do not yield reasonable 
matches (MS and RGB location, TO shape, and 
core contraction gap) to all three cluster CMDs 
due to the redder RGB. Therefore, we exclude 
the Z=0.007 solutions and only use Z=0.004 so- 
lutions for further analysis. However, it should be 
noted that NGC 2173 might be more metal-poor 
than Z=0.004 (Paper 1). For NGC 2155, OS=0.3 
isochrones produce very poor matches to the ob- 
served CMDs and thus are excluded in the further 
analysis of this study. The final isochronc solu- 
tions are plotted with the observed CMD of each 
cluster (Figure 4). Note that the shapes of all 
three isochrones for each cluster are very similar 
except for the core contraction gap. 

4- 1-2. Unresolved binary stars 

i) The effect of binaries and overshoot on 
the shape of the gap 

The effects of unresolved binary stars require 
special attention in cluster CMD studies since the 
binary sequence modifies the distribution of stars 
from the expected single star sequence. This con- 
tamination effect on the TO topology is more pro- 
found to intermediate-age clusters than to much 
younger or much older clusters since the core con- 
traction gap is prominent only in intermediate-age 
clusters. We have performed CMD synthesis in- 
cluding binary stars in order to investigate how 
they affect the distribution of stars aroimd the 
convective core contraction gap. For this test, we 
selected isochrone solutions of SL 556 with various 
BFs (0 to 30%). 

Figure 5 shows how the shape of the CMD 
around the gap is affected by binaries and OS. For 
a given BF, larger OS results in a longer MS ex- 
tending further toward the red due to the delayed 
core contraction phase. Consequently, the color 
distribution becomes much broader around the 
gap. Also, since the evolutionary time scale after 
core contraction phase is much faster with larger 
OS, the subgiant branch is less populous with in- 
creasing OS. For a given OS parameter, more stars 
are located around the gap with increasing BF. 
Binary stars make a rounded sequence, which re- 
sembles the shape of the single star sequence near 
the TO, which is » 0.75 magnitude fainter. The 
stars above the gap are most likely binaries since 



single stars cannot be located above the gap even 
with much larger OS. 

Binaries fill in the core contraction gap, which 
causes visual confusion near the TO and makes 
the derivations of age and OS parameter difficult. 
However, their sequence can still be differentiated 
from that of single stars using a detailed TO topol- 
ogy comparison. Furthermore, the number of stars 
above the core contraction gap, which are most 
likely unresolved binaries according to our analy- 
sis, can constrain the BF. 

ii) Constraining the binary fraction 
The binary frequency in a cluster CMD can be 
determined by artificial star tests using a binary 
sequence parallel to the MS (Rubenstein & Bailyn 
1997). However, this method is time-consuming 
and still has a large uncertainty in estimating the 
BF. For our clusters, which are of intermediate-age 
and expected to display a core contraction gap, the 
distribution of stars around the core contraction 
gap can be utilized for estimating the BF. 

We estimate the BF using the number ratio of 
the stars above the core contraction gap to the 
stars in the upper MS, respectively counted in two 
boxes in the CMD (Figure 6). The location of 
Box 1 is carefully selected to contain mostly bi- 
nary stars and some subgiant stars (although for 
the oldest cluster NGC 2155 the number of sub- 
giant stars increases). As mentioned in the previ- 
ous section, we have compared a synthetic CMD 
without binaries to a CMD with 30% binaries, and 
determined the best location for the binary star 
box (Box 1). The location of Box 2 is selected to 
cover the TO region including a large fraction of 
MS stars and some binaries. If the total number of 
stars is fixed, a larger value of BF naturally results 
in an increase in number of stars in Box 1, while 
the number of stars in Box 2 is not significantly 
affected by the BF. Therefore, the number ratio of 
Box 1 to Box 2 can be a good indicator of the BF 
in the observed CMD. 

The number ratio of stars in these two boxes 
for each cluster is shown in Figure 7 with a solid 
line. The uncertainties of the number ratios in the 
obscrvc^d CMDs arc estimated from the field star 
decontamination process and shown with dashed 
lines. In spite of the same box size, the observed 
number ratio of SL 556 is much greater than that 
of NGC 2173, which indicates a larger BF in 
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SL 556. This result is consistent with the pres- 
ence of a clear binary sequence parallel to the MS 
in the CMD of SL 556. In case of NGC 2155, a 
direct comparison with the other clusters is not 
possible since a Box 2 of much smaller size has 
been used due to the brighter limiting magnitude. 

We have performed 10 Monte Carlo simulations 
for each given BF and averaged the number of 
stars in each box. Dots in Figure 7 show the mean 
number ratio of the two boxes in the synthetic 
CMDs, with error bars indicating 1 a dispersion. 
Simulations show a general trend of the binary ef- 
fect on the number ratio, which increases with in- 
creasing BF as expected. According to our simula- 
tions, the number ratio varies slightly with differ- 
ent box locations and sizes, however, the trend ap- 
pears to remain the same if we keep the same boxes 
in the observed CMD and the synthetic CMD. A 
comparison between the number ratios from the 
empirical data and from the simulations suggests 
10-20% of unresolved binaries (with q > 0.7, the 
BF could be larger for any q) are required for 
all three clusters. These estimates are similar to 
those of Elson et al. (1998) for the LMC cluster 
NGC 1818 and of Rubenstein & Bailyn (1997) for 
the Galactic globular cluster NGC 6752. 

4.1.3. Overshoot from the color distribution 

around the core contraction gap 

As mentioned in § 3.1, the empirical TO mor- 
phology may constrain the amount of OS and find 
the most probable isochronc solution in Tabic 1 
and Figure 4. For a quantitative analysis, we have 
drawn Box 3 around the core contraction gap and 
analyzed the color distribution of stars in the box. 
The location of the box for each cluster, which 
contains most stars in core contraction phase, is 
shown in Figure 8. 

For all synthetic CMD models with each 
isochrone solution (OS and age), the color distri- 
butions of stars in Box 3 are plotted in Figure 9. 
Each panel contains the distribution of stars from 
the observed CMD (histogram) and from simula- 
tions with various values of BF (10%-30%), which 
are suggested in the previous section. Note that 
once scaled to the total number of stars in Box 
3, the color distribution is almost independent of 
BF, while OS affects the color distribution signif- 
icantly. In general, the distribution is broadened 
and shifted to the red with increasing OS. 



The best fits for NGC 2173 and SL 556 are mod- 
els with OS=0.2, as clearly indicated by the me- 
dian of the distribution. The various BF does not 
significantly affect the color distribution and, thus, 
OS estimation. The best OS value for NGC 2155, 
which is older than the other clusters, seems to 
lie between 0.1 and 0.2. Wo conclude that for all 
3 LMC clusters the OS parameter is around 0.2. 
The fact that OS seems to be between 0.1 and 
0.2 for NGC 2155 may indicate that OS slightly 
decreases with cluster age. A larger number of 
clusters in this age range would be necessary to 
confirm the trend. 

4.2. Mass loss and the location of the red 
clump 

The distinct feature of the red clump (RC) stars 
in the CMDs of intermediate-age clusters deserves 
considerable attention for its potential use as a 
distance indicator. However, it can be useful only 
when its characteristics (i.e., luminosity depen- 
dence on age and metallicity) are properly under- 
stood (Girardi & Salaris 2001; Grocholski & Sara- 
jedini 2002). 

An age dependence of the RC luminosity is nat- 
urally expected for the younger clusters with TO 
mass > I.8M0, for which helium ignites in a non- 
degenerate core, since the core mass of helium 
burning stars increases as the age of a cluster is 
decreasing. In contrast, in the case of stars un- 
dergoing the He-flash, the degenerate core mass is 
nearly constant and independent of age, but the 
envelope mass can play a role in controlling the RC 
luminosity since the energy output from the hy- 
drogen burning shell is decreasing with a smaller 
envelope mass. Thus, if mass loss during the RGB 
phase depends on age, there can be an age depen- 
dence in the RC luminosity of intermediate-age 
clusters. 

In this section, we derive and discuss mass loss 
for each cluster using an empirically obtained rela- 
tion between RC luminosity and RC mass. First, 
we have performed synthetic CMD simulations, 
using He-burning tracks (Yi et al. 1997), and 
measured the RC luminosity. My, as a function 
of the mean RC mass (Figure 10 top panel). It 
is shown that the median RC magnitude increases 
with larger RC mass, which is a generic feature 
of He-burning stellar tracks at this mass range 
(Vandenberg 1985; Seidel, Demarque & Weinberg 



8 



1987). Second, RC mass for each observed CMD 
is empirically derived from this luminosity-mass 
relation (points in Figure 10 top panel). The ob- 
served RC magnitude is converted to My with the 
distance modulus, and a typical error (0.1 mag.) is 
assumed considering the uncertainties in distance 
modulus. It is interesting to note that NGC 2173 
and SL 556 have a fainter RC magnitude than the 
older cluster NGC 2155, which is also shown in 
Figure 4, and consequently a smaller mean RC 
mass is inferred. Third, the amount of mass loss 
along the RGB has been estimated by subtract- 
ing the mean RC mass from the RGB tip mass 
of each isochrone. The mass loss errors are cal- 
culated from the assumed uncertainties in the RC 
magnitude. 

We note that the estimated mass loss value 
for NGC 2155 is similar to that obtained using 
mass loss efficiency, 77= 0.5 - 0.7 in Reimers' for- 
mula, which successfully reproduces the horizon- 
tal branch morphology in old globular clusters 
(Yi et al. 1999). However, the mass loss derived 
for SL 556 is too large to be compatible with 77= 
0.5-0.7 and may indicate a much higher mass loss 
efficiency in the context of Reimers' formula. To 
illustrate this, we derived RC luminosity. My, for 
clusters with various age assuming a fixed mass 
loss value, 0.2 and 0.6 Mq, respectively (Figure 10 
bottom panel). 

The origin of such a large estimate of mass loss 
in SSL 556 and disagreement with 77=0.5-0.7 re- 
mains unclear to us and may have an important 
implication on the uncertainty/variability of the 
mass loss efiiciency even in similar environments. 
The validity of mass loss estimates are subject 
to the uncertainties in distance modulus, age and 
metallicity although the color difference between 
the TO and the RGB is clear enough to indicate a 
reasonably accurate age estimate and the assumed 
metallicity provides a good match on the RGB 
slope. Further investigations are essential to draw 
more confident conclusions, especially with a bet- 
ter distance modulus constraint. 

In the case of NGC 2173, it is imclear whether 
the RC stars experienced He-flash since the RGB 
tip mass of this cluster lies near the critical mass 
at which stars develop a degenerate core, although 
it is known that the critical mass becomes smaller 
with increasing OS (Chiosi et al. 1989). If not, 
then the lower RC luminosity of NGC 2173 may 



result from the fact that RC luminosity decreases 
as RGB tip mass becomes much larger than 1.5 
Mq (see Girardi, Mermilliod & Carraro 2000). 
However, this cannot be the case for SL 556 and 
NGC 2155 since their RGB tip mass is similar or 
less than ~ 1.5 M©. 

5. Discussion: Age gap among LMC clus- 
ters 

Jensen, Mould & Reid (1988) first noticed a 
gap in the age distribution of LMC clusters, with 
none of them with ages between 4 and lOGyr. 
This age gap would be present also in metallicity, 
with younger LMC clusters having [Fe/H]~ —0.7 
and older ones having [Fe/H]2± —2.0 (Olszewski 
et al. 1991). Geisler et al. (1997) performed a 
survey of candidate old clusters and confirm the 
existence of a gap in the LMC cluster formation, 
which they situate similarly between 3 and 8 Gyr 
ago. Shortly after, Sarajedini (1998) claimed to 
have found three 4 Gyr old LMC clusters, namely 
NGC 2155, SL 556 and NGC 2121, which would 
be in the age gap. With deeper observations. Rich, 
Shara & Zurek (2001) have estimated the age of 
these clusters to be 3.2 Gyr old, and claimed that 
they are probably the first clusters to have formed 
in the younger group of the LMC clusters. Our 
age estimate of NGC 2155 is slightly younger (by 
~ 0.3 Gyr) but consistent with that of Rich et al. 
(2001), supporting the extension of the LMC age 
gap up to 3 Gyr on its intermediate-age side. 

6. Conclusion 

In this paper, we have performed detailed syn- 
thetic CMD simulations on the effect of OS and 
binary stars, compared with observed CMDs of 
three LMC clusters, NGC 2173, SL 556, and 
NGC 2155. Unresolved binary stars fill in the 
core contraction gap and make a unique sequence 
near the gap, which cannot be reproduced by sin- 
gle stars alone, even with a larger OS amount. We 
utilize the number ratios of stars around the core 
contraction gap to derive a first order estimate of 
the BF. Based on our star count analysis, all three 
clusters have a BF of 10% to 20% (with q > 0.7). 

Using the color distribution of stars aroimd the 
core contraction gap, which does not significantly 
depend on the BF, we find that the best OS pa- 
rameter for each cluster is close to 20% of the local 
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pressure scale height. 

With the best solutions of the BF and OS 
from star count and color distribution analysis, we 
present synthetic CMDs for each cluster with the 
observed CMD in Figure 11, 12, and 13, respec- 
tively. 

Prom the overall analysis, such as isochrone 

fitting, star counts, color distributions, and syn- 
thetic CMD comparison, we conclude that a mod- 
erate amount of OS, ~ 20% of the local pres- 
sure scale height, is essential to reproduce the ob- 
served shape around the core contraction gap in 
the CMDs of all three clusters, which are more 
metal-poor than the intermediate-age open clus- 
ters observed in the Galaxy, implying OS does not 
depend on metallicity at least for this metallicity 
range. With OS=0.2, our best age estimates are 
1.5, 2.1 and 2.9 Gyr for NGC 2173, SL 556 and 
NGC 2155, respectively. 

Wc constraintcd mass loss along the RGB phase 
from the observed RC luminosity of each cluster, 
utilizing the RC luminosity - mass relation de- 
rived from synthetic CMDs. We found that the 
mas loss estimate for NGC 2155 is consistent with 
the typical mass loss parameterization that works 
for the Galactic globular clusters. However, the 
SL 556 data, and perhaps the data for the other 
young cluster NGC 2173 as well, indicate a sub- 
stantially larger amount of mass loss than sug- 
gested by Reimers' formula. This may indicate 
once again the complexity of mass loss processes. 
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Fig. 1.— Isochrones with OS=0.1 and 0.3 for Z = 

0.004 and 0.007, respectively. Isochrones between 

. A o Arr^ 1 3 Gyr are plotted for two metaUicities. For each 
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macros v5 given age, isochrones with larger OS have brighter 

TOs and subgiant branches. 

Table 1: Isoclirone solutions for each cluster 



Object 


Z 


Overshoot 


Age (Gyr) 


Mrgb 


E(V-R) 


[m - M)o 


NGC 2173 


0.004 


0.1 


1.3 


1.751 


0.06 


18.5 






0.2 


1.5 


1.695 


0.06 


18.5 






0.3 


1.8 


1.630 


0.06 


18.5 


SL 556 


0.004 


0.1 


1.9 


1.556 


0.06 


18.4 






0.2 


2.1 


1.520 


0.06 


18.4 






0.3 


2.4 


1.485 


0.06 


18.4 


NGC 2155 


0.004 


0.1 


2.7 


1.396 


0.03 


18.5 






0.2 


2.9 


1.377 


0.03 


18.5 
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Fig. 2. — Synthetic CMDs before and after the 
error simulations. Left: Single stars are dis- 
tributed on top of an isochrone. Unresolved bi- 
naries are randomly chosen from single stars, and 
dispersed according to the various mass ratio, q 
(0.7 < q < 1). Right: Each star is scattered with 
a photometric error that is randomly chosen from 
the artificial stars errors table. More than 10% of 
the stars are lost due to the completeness simula- 
tion. 



Fig. 3. — A (V-R) dependence on the OS param- 
eter. We compare the 2 Gyr isochrones of 3 dif- 
ferent OS parameters, which are vertically shifted 
to clearly demonstrate the color difference. The 
color difference between the base of the RGB and 
the TO increases as a function of OS since stel- 
lar mass at the TO increases with larger OS for a 
given age. 
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Fig. 4. — Isochrone solutions for each cluster. 
Isochrones for each overshoot-age pair (see Table 
1) are plotted on top of the observed CMD. Due to 
the uncertainties in field star subtractions, much 
weight is put on the populous region of the CMD 
in the isochrone fitting process. For example, it 
is hard to define SGB and RGB on the observed 
CMD of NGC 2173 since there are only a few sub- 
giant stars (some of them are binaries and field 
stars) and many stars on the upper RGB seem to 
be field stars. The isochrone shapes of all three 
solutions are very similar except the core contrac- 
tion gap. 



Fig. 5. — The effect of overshoot and binary 
fraction. Various shapes of the core contrac- 
tion gap are presented in synthetic CMD simu- 
lations with isochrone solutions of SL 556. Sin- 
gle stars (crosses) and binaries (squares) have dis- 
tinct shapes around the gap. As OS increases, the 
MS gets redder and longer due to the delayed core 
contraction. Stars above the core contraction gap 
cannot be reproduced by single stars with even 
larger amount of OS. For a given OS parameter, 
more stars fill in the gap with increasing BF and 
make a distinctive round shape, which resembles 
the shape of single stars near the TO. 
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Fig. 6. — The location of boxes for testing binary 
fraction. Two boxes are shown in each cluster's 
CMD. From the binary simulation, the location of 
the binary star box (Box 1) is carefully selected in 
order to include most likely binaries. The lower 
box (Box 2) is selected to count MS stars with 
some binaries. We have used boxes of the same 
size, except for NGC 2155, for which the limiting 
magnitude is much brighter. 



Fig. 7. — The number ratio of Box 1 to Box 2 
for each isochrone solution. The average number 
ratio of Box 1 to Box 2 is estimated from 10 syn- 
thetic simulations with various BFs (circle). Dis- 
persion from the mean is denoted with an error 
bar. As BF increases, the number ratio increases 
since more binaries are located in Box 1 . The hori- 
zontal lines show the number ratio (solid line) from 
the observed CMD and the errors (dashed line) es- 
timated from the decontamination process. It is 
shown that 10 - 20 % of binaries (with q > 0.7) 
matches the number ratio of stars above the core 
contraction gap to stars in the upper MS. 
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Fig. 8. — The location of Box 3 near the core con- 
traction gap. A box covering most stars near the 
core contraction gap is represented in the CMD of 
each cluster. The color distribution of stars in this 
box for observed and synthetic CMDs will be an- 
alyzed to discriminate for the best OS parameter. 
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Fig. 9. — Color distribution of stars around the 
core contraction gap. The color distribution of 
the stars in Box 3 from synthetic simulations 
(lines) are compared with that of observation (his- 
togram). For each OS parameter, 3 different BF 
values arc used (BF=0: solid, BF=0.1: dotted, 
BF=0.2: dashed hno). Note that the effect of 
BF on the color distribution is insignificant. In 
contrast, larger OS tends to make a broader color 
distribution. The median color is also shifted to a 
redder color with increasing OS. Comparing with 
the observed color distribution, the best OS pa- 
rameters are determined. 
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Fig. 10.- Top: relation between the median 
RC luminosity and the mean RC mass. Median 
magnitudes of RC stars are calculated from simu- 
lated CMDs with various mean RC stellar masses 
(solid line). Using this relation, the mean RC 
mass of each cluster is determined from the ob- 
served RC magnitude. Note that the oldest clus- 
ter, NGC 2155 has the brightest RC luminosity. 
Bottom: Median RC magnitude of model predic- 
tions assuming fixed mass loss, 0.2 Mq (dashed 
line) and 0.6 Mq (dotted line) compared with that 
of observation (points). The age of all three clus- 
ters are adopted from the OS=0.2 isochrone solu- 
tion. 
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Fig, 11, Synthetic CMD of NGC 2173, Each 
synthetic CMD for a given pair of OS and BF is 
compared with the observed CMD. The synthetic 
CMD with the best solution (OS = 0.2) is empha- 
sized with a thick box. 
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Fig. 12. Synthetic CMD of SL 556. Each syn- 
thetic CMD for a given pair of OS and BF is 
compared with the observed CMD. The synthetic 
CMD with the best solution (OS = 0.2) is empha- 
sized with a thick box. 
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Fig. 13.— Synthetic CMD of NGC 2155. Each 
synthetic CMD for a given pair of OS and BF is 
compared with the observed CMD. On this plot, 
it is difficult to discriminate between OS=0.1 and 
0.2. 
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Table 2: Mass loss determined from RC luminosity 
with OS=0.2 



Object 


My 


RGB (Mq) 


RC(Me) 


ML{Mq) 


ML (??=0.5-0.7) 


NGC 2155 


0.2±0.1 


1.377 


1.262 


0.115± 0.152 


0.13-0.18 


SL 556 


0.4±0.1 


1.520 


0.937 


0.583± 0.167 


0.11-0.16 
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